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(Dated: July 9, 2001) 
The B° and B + meson lifetimes have been measured in e + e~ annihilation data collected 
in 1999 and 2000 with the BABAR detector at center-of-mass energies near the T(4S r ) resonance. 
Events are selected in which one B meson is fully reconstructed in a hadronic final state while 
the second B meson is reconstructed inclusively. A combined fit to the B° and the B + decay 
time difference distributions yields t b o = 1.546 ± 0.032 (stat) ± 0.022 (syst) ps, r s + = 1.673 ± 
0.032 (stat) ± 0.023 (syst) ps and t b+ /t b o = 1.082 ± 0.026 (stat) ± 0.012 (syst). 

PACS numbers: 13.25.Hw, 12.39.Hg 



The spectator quark model predicts that the two 
charge states of a meson with one heavy quark Q 
(Qu and Qd) have the same lifetime. Deviations from 
this simple picture are expected to be proportional 
to 1/toq |j], 0. Therefore, any lifetime differences are an- 
ticipated to be much smaller for bottom than for charm 
mesons. Various models |j], |) predict the ratio of the 
B + and B° meson [0 lifetimes to differ by up to 10% 
from unity. At present, this ratio is measured to be 
t b+/ t b° = 1-062 ± 0.029 0], with the most precise val- 
ues obtained by experiments operating near the Z and 
at hadron colliders. 

The lifetime measurements described here are based on 
a sample of approximately 23 million BB pairs recorded 
near the T{AS) resonance with the BABAR detector at 
the Stanford Linear Accelerator Center. The PEP-II 
asymmetric-energy e + e~ collider produces B + B~ and 
B°B° pairs moving along the beam axis (z direction) 
with a nominal Lorentz boost of (3j — 0.56. Hence, 
on average, the two B decay vertices are separated by 
(|Az|) = /?77g ms ct fa 270 /im, where r is cither the 
B° or B + lifetime, and 7^ ms is the Lorentz factor of the 
B mesons in the Y(AS) rest frame. This separation al- 
lows B lifetimes to be measured at the T(AS), with good 
statistical precision and systematic error sources different 
from those in previously published results. 

In this analysis, one of the B mesons in an event, de- 
noted B rec , is fully reconstructed in a variety of two-body 
charm and charmonium final states. The decay point of 
the other B in the event, B opp , is reconstructed inclu- 
sively. The decay probability distribution is given by 
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where At = t ICC — t opp is the (signed) difference of the 
proper decay times of the B mesons. The time inter- 
val At between the two B decays is determined from Az, 
including an event-by-event correction for the direction 
of the B mesons with respect to the z direction in the 
T(4S) frame. The challenge of the measurement is to dis- 
entangle the resolution in Az, 190 /im on average, from 
the effects of the B lifetime, since both contribute to the 
width of the At distribution. In the absence of back- 
ground, the measured At distribution is described by the 



probability density function (PDF) 

/+oo 
g{At'\r) K(At - At', a\d) d(At'), 
-OO 

(2) 

where 1Z is the At resolution function with parameters a, 
and a is the event-by-event error on At calculated from 
the vertex fits. An unbinncd maximum likelihood fit is 
used to extract the B° and B + lifetimes from the At dis- 
tributions for B°B° and B + B~ events. 

The BABAR detector is described in detail elsewhere [pi . 
Charged particle trajectories are measured by a combi- 
nation of a silicon vertex tracker (SVT) and a drift cham- 
ber (DCH) in a 1.5-T solenoidal field. For 1 GeV/c tracks, 
the impact parameter resolutions in z and in the trans- 
verse plane are 65 /xm and 55 /im, respectively. Photons 
and electrons are detected in the CsI(Tl) electromag- 
netic calorimeter (EMC). A ring imaging Cherenkov de- 
tector, the DIRC, is used for charged hadron identifica- 
tion. The DCH and SVT also provide ionization mea- 
surements, dE/dx, for particle identification. The in- 
strumented flux return (IFR) is segmented and contains 
resistive plate chambers to identify muons. Electron can- 
didates are required to have a ratio of EMC energy to 
track momentum, an EMC cluster shape, DCH dE/dx 
and DIRC Cherenkov angle consistent with the electron 
hypothesis. Muon candidates are required to have an 
energy deposit in the EMC consistent with the muon 
hypothesis, IFR hits located consistently on the extrapo- 
lated DCH track, and an IFR penetration in interaction 
lengths consistent with the muon hypothesis. 

B° and B + mesons are reconstructed in a sample 
of multihadron events in the modes B° — > D^~tt + , 
DW-p+, dW-o}", JfyK* and B+ -> D^°tt+, 
J/%pK + , ip(2S)K + . Multihadron events must have a 
minimum of three reconstructed charged tracks, a total 
charged and neutral energy greater than 4.5 GeV, and an 
event vertex within 0.5 cm of the beam spot [p| center 
in xy and within 6 cm in z. The event vertex is de- 
termined from all charged tracks that have an impact 
parameter with respect to the beam spot center smaller 
than 1 cm in xy and 3 cm in z. 

For 7T° candidates, pairs of photons in the EMC, each 
with more than 30MeV of energy, are selected if their 



invariant mass is within 20 MeV/c 2 of the tt° mass [0 and 
their total energy exceeds 200 MeV (100 MeV for the soft 
7T° in D* decays). A mass constraint is applied to selected 
candidates for use in the subsequent reconstruction chain. 

K® — ► 7T + 7r _ candidates are required to have an in- 
variant mass between 462 and 534MeV/c 2 . A geometri- 
cal vertex fit with \ 2 probability above 0.1% is required, 
and the transverse flight distance from the event vertex 
must be greater than 0.2 cm. 

D° candidates are reconstructed in the decay chan- 
nels K + ir~ , K + 7r~ir°, K + ir + ir~TT~ and K^tt + tt~ and 
D~ candidates in the decay channels K + it~it~ and 
K®tt~. Kaons from D~ decays and charged daughters 
from D° — > K + tt~ are required to have momenta greater 
than 200MeV/c. All other charged D daughters are re- 
quired to have momenta greater than 150MeV/c. For 



is required to be greater than 0.1%. Positive identifi- 
cation of kaons is required for modes with higher back- 



D° -> K+i 



we only reconstruct the dominant reso- 



nant mode, D° — * K + p~ , followed by p~ 
Tr~ir° mass is required to lie within 150MeV/c 2 of the 
p mass H and the angle between the ir~ and D° in the 
p rest frame, 6* D0 , must satisfy |cos0£o„.| > 0.4. All 
D° and D~ candidates are required to have a momentum 
greater than 1.3GeV/c in the T(4S) frame, an invariant 
mass within 3a of the nominal value H and a geomet- 
rical vertex fit with a \ 2 probability greater than 0.1%. 
A mass constraint is applied to selected D candidates. 

Charged and neutral D* candidates are formed by 
combining a D° with a 7r~ or 7r°. The momentum of 
the pion in the T(4S) frame is required to be less than 
450 MeV/c. The soft 7r~ is constrained to originate from 
the beam spot when the D*~ vertex is fit. After the 
mass constraint to the D° daughter, D* candidates with 
m(D°Tr) within 2.5a of the nominal mass M for D*~ , or 
within 4(7 of the nominal mass Gj for D*°, are selected. 

Candidates for leptonic decays of charmonium mesons 
must have at least one decay product positively identified 
as an electron or a muon. If it traverses the calorimeter, 
the second muon must be consistent with a minimum 
ionizing particle. J/ip candidates are required to lie in 
the invariant mass interval 2.95 (3.06) to 3.14GeV/c 2 for 
the e + e~ (p + p~) channel. The e + e~ (p + p~) invariant 
mass of ip(2S) candidates must be between 3.44 (3.64) 
and 3.74 GeV/c 2 . A mass constraint is applied to selected 
candidates. ip(2S) — ► J/07r + 7r~ candidates are selected 
if the 7r + 7r~ mass is between 0.4 and 0.6 GeV/c 2 and the 
ijj(2S) mass is within 15MeV/c 2 of the nominal value |jj. 
All ip(2S) candidates must have momenta between 1.0 
and 1.6GeV/c in the T(4S) rest frame. 

B candidates are formed by combining a D^*' , J/ip or 



ip{2S) candidate with a ir~ 



(at 



K*° (K* Q -> K+tt-) or K+ candidate that has a mo- 
mentum larger than 500 MeV/c in the T(4S) frame. For 
B° — » D^*'~ p + , the 7T° from the p + decay must have an 
energy greater than 300 MeV. For B° -> DW'aJ, the af 
must have an invariant mass between 1.0 and 1.6 GeV/c 2 , 
and the \ 2 probability of a vertex fit of the af candidate 
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with D° 
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ground, such as B + 

Continuum e + e~ — ► qq background is rejected by re- 
quiring the normalized second Fox- Wolfram moment o 
for the event to be less than 0.5. Further suppres- 
sion is achieved by a mode-dependent restriction on the 
angle between the B Tec and i? pp thrust axes in the 
T(45) frame. 

B° and B + candidates are identified on the basis 
of the difference AE between the reconstructed energy 
and the beam energy t/s/2 in the T(4S) frame, and 
the beam-energy substituted mass tties calculated from 
y/s/2 and the reconstructed momentum of the candidate. 
B candidates are selected with toes > 5.2 GeV/c 2 and 
|AE| < 3ct A e, where ct A e (10 to 30 MeV) is the mea- 
sured resolution for each decay mode. 

The decay position of the B rcc candidate is determined 
by requiring convergence of a vertex fit, where in addition 
the masses of the D mesons are constrained to their nom- 
inal values m. Precisions between 60 and 100 pm rms for 
the B rcc decay position in z and in the transverse plane 
are achieved, depending on the decay mode. 

The vertex of the -B pp is determined from all tracks 
in the event after removing those associated with the 
B rcc candidate. Tracks from photon conversion candi- 
dates are rejected. Daughter tracks from K° s or A can- 
didates are replaced by the neutral parents. An addi- 
tional constraint is imposed on the -B pp vertex using 
the By- cc vertex and three-momentum, the beam spot po- 
sition, and the average T(4S) momentum. To reduce 
the bias in the forward z direction from charm decay 
tracks, the track with the largest contribution to the ver- 
tex x 2 , if above 6, is removed and the fit iterated until 
no track fails this requirement. Events are required to 
have at least 2 tracks remaining in the B opp vertex, an 
error on Az smaller than 400 pm and |Az| < 3000 pm. 
The precision achieved on Az, 190 pm rms on average, is 
dominated by the resolution on the B opp vertex. A re- 
maining bias of —35 pm due to charm decays on the 
B opp side is observed. We require |Ai| < 18 ps and find 
6064± 70 B° and 6336± 63 B+ signal events in a ±2.5cr 
{a = 2.7 MeV/ c 2 and 2.6MeV/c 2 , respectively) window 
around the toes peak above a small background (~ 10%). 
The toes distributions for the final samples are shown in 
Fig. along with the results of a fit with a Gaussian 
distribution for the signal and an ARGUS background 
function M. 

As already noted, the modeling of the resolution func- 
tion 1Z is a crucial element of the B lifetime measure- 
ments. Studies with both Monte Carlo simulation and 
data show that the sum of a zero-mean Gaussian dis- 
tribution and its convolution with a decay exponential 
provides a good trade-off between different sources of un- 
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FIG. 1: tubs distributions of the selected neutral (top) and 
charged (bottom) _B rcc candidates. 
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where St is the difference between the measured and true 
At values. The model parameters a are the fraction h 
in the core Gaussian component, a scale factor s for the 
per-event errors a, and the factor k, in the effective time 
constant ko~ of the exponential that accounts for the ef- 
fect of charm decays. Monte Carlo studies show that 
the parameters a obtained for different decay modes are 
compatible, as expected for At resolution dominated by 
the -B pp vertex. The resolution functions for B° and B + 
mesons differ slightly because the B opp decays to a differ- 
ent admixture of D~ and D° mesons. The difference is 
not significant given the present data sample size. Hence 
a single set of resolution function parameters is used for 
both B° and B + in the lifetime fits, and a small correc- 
tion discussed later is applied to the results. While the 
resolution function 1Z describes almost all events, incor- 
rectly measured outlier events are modeled separately as 
discussed below. 

The unbinned maximum likelihood fit for the B life- 
times uses all events with to-es > 5.2GeV/c 2 . The 
probability p^ 9 for event i to be signal with At dis- 
tribution Q, defined in Eq. |2j, is estimated from the 
toes fit (Fig. |l|) and the toes value of the B T0C candidate. 
Each event i then samples a PDF that includes signal, 
background, and outlier components: 
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The background At distribution, B, for each B species is 
modeled by the sum of a prompt component and a life- 
time component convoluted with a resolution function of 
the form given in Eq. [|, but with a separate set of pa- 
rameters. The fraction of non-prompt background, its 
effective lifetime and the background resolution param- 
eters are determined separately for charged and neutral 
B mesons. Signal and background outlier events have an 
assumed At behavior O given by a Gaussian distribution 
with zero mean and a fixed 10 ps width. The fractions 
of outliers in signal and background are determined sep- 
arately in the lifetime fit. 

Since the same resolution function is used for neutral 
and charged B mesons, the fitting procedure maximizes 
the log-likelihood function ln£ formed from the sum of 
two terms, one for each B meson species, with common 
parameters a for 1Z: 

]n£ = £ HT{At l+ ,a %+ ,ptX\r B+ ; a, b+,fo7i + Jou?' + )\ 
+ ^]n[^(Atio,aio,p s ^ 9 \T B o-,a,bo,fout Jott' )]- ( 5 ) 

iO 

The likelihood fit involves 19 free parameters. The pa- 
rameter t b + is replaced with t b + — r ■ t b o to estimate 
the statistical error on the lifetime ratio r. The lifetime 
values were kept hidden until the event selection and the 
At reconstruction method, as well as the fitting proce- 
dure, were finalized and the systematic errors were de- 
termined. 

The fit results, after small corrections discussed below, 
are t b o = 1.546 ± 0.032 ps, t b + = 1.673 ± 0.032 ps and 
t b +/t b o = 1.082 ± 0.026, where the errors are statisti- 
cal only. The resolution parameters a (h = 0.69 ± 0.07, 
s = 1.21 ± 0.07 and n = 1.04 ± 0.24) are consistent with 
those found in a Monte Carlo simulation that includes 
detector alignment effects. The fitted outlier fractions 
in the B+ and B° signals are both 0.2l°;|%. Figure | 
shows the results of the fit superimposed on the observed 
At distributions for B° and B + events within 2.5 stan- 
dard deviations of the B mass in Toes- 

Table |l| summarizes the systematic uncertainties on 
the lifetime results. The full analysis chain, including 
event reconstruction and selection, has been tested with 
Monte Carlo simulation. The statistical precision on the 
consistency between the generated and fitted lifetimes is 
assigned as a systematic error. The resolution param- 
eters a are determined from the data by the fit, con- 
tributing ±0.017 ps in quadrature to the statistical error 
of the individual lifetime results. Thus, a large part of 
the At resolution uncertainty is included in the statistical 
error. Residual systematic uncertainties are attributed to 
limited flexibility of the resolution model. These contri- 
butions have been estimated by comparing results with 
different parametrizations. We correct our measurements 
for the small positive (negative) bias on the B° (B + ) life- 
time due to differences in the At resolution functions for 
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FIG. 2: At distribution for the B° (top) and B + (bottom) 
events within 2.5a of the B mass in tties- The results of 
the fit are superimposed on the data. The single-hatched 
areas are the background components B and the cross-hatched 
areas represent the outlier contributions. The probability of 
obtaining a lower likelihood, evaluated with a Monte Carlo 
technique, is 7.3%. 



TABLE I: Summary of the systematic uncertainties. 



Effect 


8(t b o) 


S(t b+ ) 


S(t b + /t b o ) 




(PS) 


(PS) 




MC statistics 


0.009 


0.007 


0.006 


TL parametrization 


0.008 


0.004 


0.003 


same 1Z for B° and B + 


0.004 


0.005 


0.006 


Beam spot, ps roc 


0.002 


0.002 


cancels 


At outliers 


0.011 


0.011 


0.005 


SVT alignment 


0.008 


0.008 


cancels 


z scale 


0.008 


0.008 


cancels 


Az to At conversion 


0.006 


0.006 


cancels 


Signal probability 


0.003 


0.003 


0.003 


Background modeling 


0.005 


0.011 


0.005 


Total in quadrature 


0.022 


0.023 


0.012 



B and B + mesons arising from their decays to a differ- 
ent admixture of D~ and D° mesons and estimated with 
a high-statistics Monte Carlo sample. The size of the 
correction is assigned as a systematic error. A small sys- 
tematic error results from uncertainties on the beam spot 

position and vertical si?p and the R rc - momentum ver- 



experiments with sets of fixed values for the outlier width 
and mean were simulated and subjected to the nominal 
lifetime fit. The largest observed bias is taken as system- 
atic uncertainty. Additional systematic uncertainties are 
due to the SVT alignment. The z length scale was deter- 
mined to better than 0.5% from secondary interactions 
in a beam pipe section of known length. Approximations 
in the calculation of At from Az and the uncertainty 
on the boost lead to small systematic errors. The er- 
rors on the toes fit parameters are used to determine 
the uncertainty on p SZ9 and the corresponding system- 
atic error. The main systematic uncertainties related to 
backgrounds arise from changes in the background com- 
position as a function of toes • An additional contribution 
arises from a 1-2% B° contamination of the B + signal 
sample and vice versa. We use Monte Carlo simulation 
to correct for these background effects and assign the 
sum in quadrature of the corrections as systematic un- 
certainty. 

In summary, the B° and B + meson lifetimes and their 
ratio have been determined to be: 

r B a = 1.546 ± 0.032 (stat) ± 0.022 (syst) ps, 
t b + = 1.673 ± 0.032 (stat) ± 0.023 (syst) ps, and 
t b +/t b o = 1.082 ±0.026 (stat) ±0.012 (syst). 

These are the most precise measurements to date, and 
they are consistent with the current world averages. 
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